Spin texture on the Fermi surface of tensile strained HgTe 
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We present ab initio and k • p calculations of the spin texture on the Fermi surface of tensile 
strained HgTe, which is obtained by stretching the zincblende lattice along the (111) axis. Tensile 
strained HgTe is a semimetal with pointlike accidental degeneracies between a mirror symmetry 
protected twofold degenerate band and two nondegenerate bands near the Fermi level. The Fermi 
surface consists of two ellipsoids which contact at the point where the Fermi level crosses the twofold 
degenerate band along the (111) axis. However, the spin texture of occupied states indicates that 
neither ellipsoid carries a compensating Chern number. Consequently, the spin texture is locked 
in the plane perpendicular to the (111) axis, exhibits a nonzero winding number in that plane, 
and changes winding number from one end of the Fermi ellipsoids to the other. The change in the 
winding of the spin texture suggests the existence of singular points. An ordered alloy of HgTe with 
ZnTe has the same effect as stretching the zincblende lattice in the (111) direction. We present ab 
initio calculations of ordered Hg x Zni_ x Te that confirm the existence of a spin texture locked in a 
2D plane on the Fermi surface with different winding numbers on either end. 



Mercury telluride (HgTe) is a zero band-gap zincblende 
semiconductor with an inverted ordering of states at 
where the band gap and inversion can be controlled by 
alloying with CdTe. HgTe and CdTe are both zincblende 
materials, with the crucial difference that the stronger 
spin-orbit coupling of Hg causes the states at V to in- 
vert relative to each other [Fig. [I] . Therefore HgTe/CdTe 
quantum wells exhibit the quantum spin Hall effeclP. The 
electronic properties of alloys of HgTe have been stud- 
ied extensively for many decades. Recently there has 
been a surge of interest in so-called topological (Weyl) 
scmimetals: three dimensional materials which exhibit 
pointlike degeneracies (Weyl points^) between bulk con- 
duction and valence bands and exhibit topologically pro- 
tected surface modes^^. These topological Weyl points 
are low symmetry descendants of a so-called Dirac point, 
which is described by a massless four-band Dirac-like 
HamiltoniarP. The unique feature of a Dirac semimetal 
is a Fermi surface that surrounds discrete Dirac points, 
around which all bands disperse linearly in all directions. 
It turns out that HgTe is almost a Dirac semimetal: at 
r two conduction states are degenerate with two valence 
states; all four bands disperse linearly along the (110) and 
(100) axes; however, mirror and time reversal symmetry 
guarantee that along the (111) axis, two of the four bands 



must disperse quadratically [Refs. loTH and Fig. 2 (a 



Figure 2(a) shows another interesting feature of the 
bands at T. In addition to the symmetry guaranteed four- 
fold degeneracy, there are additional pointlike accidental 
degeneracies along the (111) axis between two conduc- 
tion states and one valence state. The rest of the band 
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structure is completely gapped, so the Fermi surface is 
confined to a small region around T [Fig. [2(b)] . It is nat- 
ural to ask what perturbations could turn HgTe into a 
Dirac semimetal. Since the quadratic dispersion of bands 
near T is enforced by mirror and time reversal symmetry, 
any perturbation must break either one of these symme- 
tries. A time reversal symmetry breaking Zeeman field 
breaks at least one rotation symmetry of the tetrahedral 
group splitting the fourfold degeneracy at T so the system 
either gaps or develops additional Fermi surface. Break- 
ing mirror symmetry maintains the fourfold degeneracy 
at r, but the Dirac point develops a nonzero Chern num- 
ber due to the threefold rotation symmetry of zincblende. 
This inevitably leads to additional Fermi surface as well. 
Hence, zincblende materials cannot be Dirac semimet- 
al:P. A third option is to break, for instance, a fourfold 
rotation symmetry by straining the lattice. Such a pro- 
cedure will split the fourfold degeneracy at T and will 
either gap the band structure completely or else shift the 
accidental degeneracy away from T along the (111) axis. 



Reference [H] predicted that compressive strain along 
the (111) axis will gap HgTe into a topological insulator, 
since the states at T are already inverted [Fig. 3(a)] . The 
inverted ordering of bands at T is a direct consequence of 
spin-orbit coupling. To understand how it comes about, 
we compare the ordering of states in HgTe with and with- 
out spin [Fig.[l]. At the Fermi level in spinless HgTe, the 
valence states span a p-type representation T4 of space 
group 216. Slightly above in energy, there exists an 
s-type band which belongs to the symmetric represen- 
tation Ti. When spin-orbit coupling is introduced, the 
representation T4 splits into a £>3/2-type representation 
r% and a Pi/2-type representation F7, whereas Ti turns 
into an Si/2-type representation Tq. Since Hg is a heavy 
atom, its contribution to the spin-orbit coupling is strong 
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FIG. 1: Schematic illustration of band inversion in HgTe due 
to strong spin-orbit coupling. Twofold degenerate bands are 
drawn as thick lines, whereas nondegenerate bands are drawn 
as thin lines. Left: Without spin, HgTe/CdTe have an s-type 
conduction band and a p-type valence band. Middle: The 
spin-orbit coupling in HgTe is strong enough to push the Si/ 2 - 
type representation T@ below the P3/2-type representation T%. 
Right: Normal ordering of states in CdTe, which has a weaker 
spin-orbit coupling as compared with HgTe. The S!/ 2 -type 
representation Fs appears above the P3/ 2 -type representation 
Tg, as expected from the spinless case. 



enough to push the Si/2-type band T 6 below r 8 , resulting 
in an inverted band structure. This is in contrast with 
a normal ordered zincblende semiconductor (e.g. CdTe), 
where Tq appears above r 8 in energy. Reference[2]demon- 
strated that if T 8 is split and the inversion of T 6 relative 
to the split band r 8 is controlled by a "mass" parameter, 
the insulating phases that ensue have a non-trivial Z2 
index +1 in the inverted regime, versus a trivial index 
in the normal regime. The splitting of r 8 was achieved 
in a HgTe/CdTe quantum well structure, while the tun- 
ing of Tq above or below the split representation T 8 was 
controlled by the well thickness, leading to the first ex- 
perimental observation of the quantum spin Hall effect 9 . 

In this paper, we show that whereas compressive strain 
along the (111) axis gaps the system entirely [Fig. |^a)] , 
tensile strain along the same axis shifts the accidental 
degeneracies further along the (111) axis after splitting 
^[b)| . Tensile strained HgTe remains a semimetal 



r 8 [Fi g _ 

where the dispersion of bands around the accidental de- 
generacies is linear in all directions [Fig. [4]. But there 
is a caveat. Each accidental degeneracy involves three 
states so it is neither a Dirac point nor a Weyl point. 
As discussed earlier, this is due to the presence of mirror 
symmetry about a plane that contains the (111) axis. If 
we could break the mirror symmetry, the twofold degen- 
erate band in Fig. |2(a) would split and the system will 
develop two Weyl points [Fig. [5] . 

A Weyl point carries a nonzero Chern number (charge) 
±1 and cannot be annihilated unless brought in contact 
with another Weyl point of opposite charge. In other 
words, the spin texture of valence states around a Weyl 
point has nonzero divergence; if interpreted as a mag- 
netic field, the spin texture would correspond to a mag- 



netic monopole with charge ±1. However mirror sym- 
metry constrains all points on its mirror plane to have 
zero Chern number. Hence Weyl points cannot exist on 
mirror planes, and so the integral of the normal compo- 
nent of spins on a closed surface that encloses points on 
a mirror plane must be zero. 

This has important implications for the Fermi surface 
of tensile strained HgTe [Fig. |4j. The two Fermi ellip- 
soids contact at a point due to mirror symmetry. They 
can be separated only if mirror symmetry is broken. Con- 
sequently their corresponding spin textures will develop 
a nonzero divergence. We ask what is the nature of 
the spin texture on the ellipsoidal Fermi surface before 
mirror symmetry is broken. It turns out that spin tex- 
ture on both ellipsoids is locked in a plane perpendicular 
to the (111) axis and exhibits a nonzero winding num- 
ber about that axis. Furthermore the winding number 
changes as one sweeps across the Fermi ellipsoids from 
one end to the other and inevitably exhibits spin singu- 
larities on the Fermi surface where the winding number 
must change [Fig. [6] . 

We alloy HgTe with ZnTe by replacing half the Hg 
atoms with Zn atoms. The size mismatch between Hg 
and Zn causes the lattice to relax so that the Zn and Hg 
atoms move off center towards the Te atoms. This has 
the same effect as tensile strain along the (111) axis in 
zincblende [Fig. ff\. Indeed, ab initio calculations indicate 
that the spin texture on the Fermi surface of Hg x Zn!_ x Te 
(x = 0.5) is locked in a two dimensional plane and ex- 
hibits spin singular points identical to tensile strained 
HgTe [Fig. [I]. However it is likely that Density Func- 
tional Theory (DFT) overestimatealSlthe amount of band 
inversion for x = 0.5 given that Hg x Cdi_ x Te is a band 
insulator for x < 0.8 11 . Nonetheless, for small enough 
concentrations of Zn, this symmetry protected spin tex- 
ture may indeed by realized physically. It is difficult to 
span the full range of intermediate compositions within 
DFT to make a prediction of the composition where the 
inversion occurs. 

In Section |TJ we describe the Fermi surface of ten- 
sile strained HgTe. In Section [TTJ we derive a simple 
k • p model to understand band dispersion in tensile 
strained HgTe. In Section |III| we derive the spin tex- 
ture at the k • p level and provide ab initio calculations 
of Hgo.5Zno.5Te in support of the predictions of the k • p 
theory. Finally we conclude with some brief remarks. 



I. FERMI SURFACE 

The band structure of HgTe has been of considerable 
interest in recent years because HgTe alloyed with CdTe 
exhibits a tunable direct band-gap from to 1.5 eV at 
r. We focus our attention to the four dimensional irre- 
ducible representation T 8 at the Fermi level. The point 
group elements of the little group at T belong to T4 
(the tetrahedral group with mirror planes). Since the 
zincblende lattice lacks inversion symmetry, the matrix 
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FIG. 2: Band structure of HgTe. (a) Energy bands near the Fermi level as a function of momentum along the (111) axis. At V 
the four degenerate states span the irreducible representation Tg of the group Td- are irreducible representations of the 

group Csv . The twofold degenerate (solid) band labeled Ag disperses quadratically to lowest order, whereas the nondegenerate 
bands A4 and A5 disperse linearly. There are threefold accidental degeneracies at k = ±0.004(1, 1, 1) because the bands A4 and 
A5 have negative mass, whereas Ag has positive mass, (b) Band structure of HgTe as a function of crystal momentum across 



the entire Brillouin zone. The valence and conduction bands are gapped everywhere except near F 
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FIG. 3: I (a) I Band structure of HgTe under compressive strain along the (111) axis. The fourfold degenerate representation 
Fg is split, and the material becomes a topological insulator^, (b) Band structure of HgTe under tensile strain along the (111) 
axis. Although Fg splits under tensile strain as well, the accidental degeneracy along the (111) axis shifts away from F and the 
material remains semimetallic. 



elements (i/>i\p\il>j) , where span Tg, are non-zero. 

Therefore, the four bands degenerate at V can, in prin- 
ciple, disperse linearly in all directions around F^. How- 
ever, along the (111) axis, T% splits into a two dimensional 
representation Aq and two one dimensional representa- 
tions A4 and A 5 of the little group along that axis. Since 
r is invariant under time reversal symmetry, A 6 is con- 
strained to be flat, whereas A4 and A 5 must disperse with 
non-zero slope. Along other directions such as (110) and 
(100), Tg splits into either four singlets or two doublets, 



and none of the bands are constrained to be fiat to linear 
order in iPEl. This is why the four band model describ- 
ing states {|r 8 , 3/2), |r 8 , -3/2), |r 6 , 1/2), |r 6) -1/2)} in 
Ref. [5] was written as a Dirac Hamiltonian in two dimen- 
sions. 

It turns out that the band s belonging to A4 and A5 
have negative mass [Fig. 2(b) . Therefore, to C(k 2 ) and 



above, the bands labeled A4 and A5 must bend back be- 
low the Fermi level and cross A 6 somewhere along the 
(111) axis. High resolution ab initio calculations confirm 
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that these degeneracies indeed occur very close to T at 
± 0.004(1, 1, 1) [Fig. 2(a)| . Since this occurs so close to 



r, the Fermi surface of HgTe is only effectively pointlikc. 

Under tensile strain along the (111) axis, the represen- 
tation Tg splits, because the tetrahedral group T<j reduces 



to the trigonal group C 3v = {1, 



, y, xy, x 2 y}. Here x 



corresponds to threefold rotation symmetry about the 
(111) axis while {y, xy, x 2 y} represent the three mirror 
planes that contain the (111) axis and are rotated by 
2tt/3 with respect to each other. Since tensile strain pre- 
serves the little group along the line from T to L, the 
representation Aq does not split. Therefore, the acciden- 
tal degeneracies between A4/A5 and A@ only shift further 
away from T, and HgTe remains a semimetal. 

Figure |4] shows the band structure along the (111) di- 
rection for 8% tensile strain in HgTe. The set of mo- 
menta corresponding to a fixed energy form two ellip- 
soidal surfaces that contact at a point. This implies that 
the bands A^As and A§ disperse linearly in all direc- 
tions around the accidental degeneracies. Therefore, un- 
der tensile strain the Fermi surface of HgTe grows from 
a point to a set of four ellipsoids, in two sets of Kramers 
paired momenta in the Brillouin zone. The two ellipsoids 
on one side of T are connected to each other at the point 
k along the (111) axis where the Fermi level crosses the 
band A6. 

We note that 8% tensile strain corresponds to 4% 
strain in-plane, which is relatively large. This effect en- 
sues as soon as a perturbation breaking the fourfold sym- 
metry of zincblcnde turns on; there is no critical value 
of strain below which the unusual spin texture of the 
Fermi surface would disappear. While 8% tensile strain, 
which corresponds to 4% strain in-plane, is quite high, 
the ellipsoidal fermi surfaces ought to be observable at 
reduced strains; at 2% tensile strain the energy difference 
between the threefold degeneracies is still 18meV, down 
from 25meV. Even so, very high strains may be achiev- 
able through epitaxial growth of HgTe on an appropriate 
substrate^. 

Alternatively, strain can be introduced chemically. Al- 
loying HgTe with ZnTe introduces a size mismatch be- 
tween Hg and Zn, causing both Zn and Hg atoms to move 
off center toward a Te atom; if all the Zn atoms move in 
the same direction, the symmetry breaking is the same 
as that created by tensile strain. A similar effect has also 
been reported in Cd x Zni_ x Tc 13 , which shows rhombohc- 
dral distortions and Zn off-centering for a wide range of 
Zn concentration. Figure [7] shows the band structure of 
ordered Hgo.5Zno.5Te, which is qualitatively the same as 
that of strained HgTe. The difference lies in the shape 
and size of the Fermi surface and the Fermi velocities in 
the two materials. As a caveat, it is important to note 
that DFT overestimates the magnitude of band inver- 
siorPH; Hg x Cdi_ x Te for x < 0.8 is a band insulato^, 
and Hgo.5Zno.5Te likely to be gapped in actuality. How- 
ever, calculations for Hgo.75Zno.25Te suggest that signif- 
icant symmetry-breaking distortion persists for lower Zn 
concentration, consistent with the observations of similar 




FIG. 4: Band dispersion in tensile strained HgTe at the Fermi 
level as a function of momentum along the (111) axis. The 
twofold degenerate band (drawn as a thick line) is acciden- 
tally degenerate with both the nondegenerate bands (drawn 
as thin lines) within a very small energy range. Band disper- 
sion around the accidental degeneracies is linear to leading 
order in the (111) direction. Right: Schematic illustration of 
the Fermi surface. The Fermi surface for various choices of 
Fermi energy (drawn in dashed lines) near the accidental de- 
generacies consists of two ellipsoids which contact at a point. 
This indicates that the bands disperse linearly in transverse 
directions as well around the accidental degeneracies. 



behavior in Cd x Zni_ x Te^. Therefore, if some ordered 
compound Hg x Zni_ x Te exists with x such that distor- 
tion due to Zn substitution sufficiently splits the degen- 
eracy at the Fermi level while still allowing band inver- 
sion, then it would possess an experimentally observable 
ellipsoidal Fermi surface. We assert that a time rever- 
sal and mirror symmetry breaking Zeeman field would 
split the degenerate band Aq and turn Hg x Zni_ x Te into 
a Weyl semimetal. This is confirmed by our k • p the- 
ory in Section [H] A Weyl semimetal obtained by a Zee- 
man perturbation to Hg x Zni_ x Te will have interesting 
phcnomcnological consequences such as Fermi-arc surface 
states^ and pressure induced anomalous Hall effeclP. 

While some experimental studies have been performed 
on Hgi-^Zn^Te, it is difficult to draw any direct conclu- 
sions without more detailed structural characterization 
of the sampleJSESl. We performed a phonon analysis 
of tensile strained HgTe and Hgo.5Zno.5Te and found no 
structural instability. 



II. K P THEORY 

We model the low energy physics of the states that 
span the representations A4, A5 and Ag along the (111) 
axis in terms of a four band k • p theory around 
the accidental degeneracies. The states which belong 
to A 6 are \p, ±|) whereas the ones belonging to A4 
and A 5 are \p, ±5). The symmetry group is C% v — 
{l,x,x 2 ,y,xy,x 2 y} with the group multiplication rules 
x 3 = y 2 = l,yx = x 2 y in addition to time reversal sym- 
metry. Here x represents a threefold rotation symmetry 
while y represents a mirror symmetry. The accidental de- 
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FIG. 5: (a) Energy bands plotted as a function of momentum 
in the k z direction in the kp theory. There is a twofold degen- 
erate band (thick line) intersecting two nondegenerate bands 
(thin lines) analogous to Fig. [4] (b) Mirror and time rever- 
sal symmetry breaking perturbation proportional to J z splits 
the twofold degenerate band into two nondegenerate bands, 
giving rise to four Weyl points, (c) Rotation and mirror sym- 
metry breaking perturbation proportional to t z g x gaps the 
system entirely. 



generacies occur at points ko 7^ along the line from T 
to L, which are not time reversal symmetric. So the time 
reversal operator maps the k • p Hamiltonian around 
ko to its time reversed partner around — ko. 

We choose an angular momentum basis for the k • p 
theory. Since we are dealing with four half-integer spin 
states, we require a double valued representation of C$ v , 
such that x and y take the form of matrices R and 
M belonging to SU(4) that satisfy R 3 = M 2 = -1. 
These operators satisfy the group multiplication rules 
and are consistent with the interpretation of the oper- 
ators M,RM,R M as mirror symmetries (all of which 
square to —1) for particles with half integer spin. 

Orienting the k z direction along the (111) axis, the 
threefold rotation operator is diagonal and can be written 
as R = exp(i J z 27r/3) where J z is the angular momentum 
operator. The mirror operator in this basis takes the form 
M = iT y (£) d x . f and a denote Pauli matrices. Tensor 
products of Pauli matrices can be used as basis vectors 
for 4x4 hermitian matrices. Neither f nor a individually 
corresponds to a real or pseudo spin degree of freedom. 
In what follows, we will write r, ® <7j as Ti<jj where the 
tensor product is implied. In the absence of inversion 
symmetry V, the local antiunitary operator VO is not a 
symmetry. VQ is a local operator because both V and 
send k to — k. In the angular momentum basis, it can 
be written as V& = ir x a y K, where K is the complex 
conjugation operator. 

A Hamiltonian linear in k which respects these sym- 
metries is, 



where 



ff(k) = k x r x + k y T y + k z T z <T z + k (1) 
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is chosen to break the local antiunitary symmetry "P0. 
With this choice of k, time reversal symmetry is pre- 
served while inversion is broken if we require that 



QH(k)S^ 1 — r x k x + Tyky + r z a z k z + k t is the effective 
k • p Hamiltonian around the time reversed location of 
the accidental degeneracies. Such a constraint is allowed 
because the k • p theory is localized around a Brillouin 
zone momentum ko 7^ along the (111) axis. 

Figure[5] shows the energy spectrum of the Hamiltonian 
in Eq. (fTT along the k z axis. The degenerate band Aq 
splits to give four Weyl points under a mirror and time 
reversal symmetry breaking Zeeman term proportional to 
a z or 3 z where 3 z is the angular momentum operator of 
p 3 / 2 states. Fig.[5|b) shows the splitting of Aq to give four 
Weyl points due to a perturbation proportional to 3 Z . 
A rotation and mirror symmetry breaking perturbation 
such as T z a x gaps the system entirely [Fig. [jjjc)]. 

It is possible to split Ag into two parallel bands as in 
Fig. j5^b) by a mirror preserving term r y a x which breaks 
time reversal and rotation symmetry. This however does 
not lead to Weyl points; even though the degeneracies 
involve two states, the bands develop a quadratic disper- 
sion in one direction. Mirror symmetry about a plane 
constrains all points k on the plane to have zero Chern 
number. This is explained as follows: The Chern number 
of a certain point k in the Brillouin zone is given by an in- 
teger n = l/(2ni) <f s TrF-n where Fij = Vk X (tpi\ VklV'j) 
is the Berry curvature, h is a unit vector normal to the 
surface S that contains k, and the trace includes only 
valence states \tpi), i = 1, ..N. A mirror symmetry about 
a plane that contains k reverses the orientation of the 
surface S while the integrand is the same for opposing 
points on S. Hence the Chern number of all points on a 
mirror plane must be zero. Therefore, as long as mirror 
symmetry is present, tensile strained HgTe cannot host 
Weyl points along the (111) axis. 



III. SPIN TEXTURE 



Since the dispersion of bands around the accidental de- 
generacies in Figs. [4] and 5ja) is linear in k, it is natural 
to ask if the accidental degeneracies have an associated 
Chern number. However since the (111) axis belongs to a 
mirror plane, the Chern number of all points on that axis 
should be zero. In other words, the surface integral of the 
normal component of the spin texture over a closed sur- 
face enclosing the accidental degeneracies will be zero. 
The ellipsoidal Fermi surface in Hgo.5Zno.5Te also en- 
closes the accidental degeneracies, as shown in Fig. [4] 
At the k • p level, this corresponds to two spheres in k 
space parametrized as, Si : £ 2 = k 2 + k y + (l — k z ) 2 and 
S 2 ■ {£ — l) 2 = k\ + k 2 + k 2 where £ is the Fermi energy. 

The average spin field is given by {ipi\S\ipi) where \ipi) is a 
valence state on one of the Fermi ellipsoids parametrized 
by £ and S is the spin operator for the p 3 / 2 states. In 
the k • p theory, this evaluates to, 
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(S)* = a i (£ ) k z ) [sin(0)i + cos 

+ [sin (20) i- cos (20) j 



(2) 



where the subscript i = 1,2 identifies one of the Fermi 
ellipsoids Sj and is the azimuthal angle in the k x k y 
plane that parametrizes the ellipsoids. The coefficients 



are, 



ai(8, k z ) 
a 2 (£, k z ) 



U(£ + l-k z ) 
V 3(£ - 1 + k z ) 

U(£-l-k z ) 
3(£-l + k z ) 



,Pi(£,k z 



,p2(£,k z ) 



ol\{£, k z 



a 2 (£,k z ) 2 



The three mirror planes of the group C^ v , all of which 
contain the k z axis, lock the spin expectation values in 
the k x k y plane. Figure [6] shows how the spin texture 
evolves as a function of k z around the Fermi ellipsoids. 
For fixed £, k z ranges between 1 — \£\ and 1 + \£ | on Si 
and between — 11 — £ | and |1 — £\ on S 2 . When k z is close 
to 1— £ on either ellipsoid, the term with the coefficient ft 
dominates and the spin winding number is —2, whereas 
when k z is close to £ + 1 on Si or £ — 1 on S 2 the term 
with the coefficent ccj dominates, and the spin winding 
number is +1. Therefore, the spin winding number must 
change abruptly from one end of the Fermi ellipsoids to 
the other, since the spins cannot develop a component in 
the k z direction due to mirror symmetry. Hence, there 
must exist singular points on the Fermi ellipsoids where 
the spin expectation values vanish. Figure [6] shows the 
points on the Fermi ellipsoids where the spin expectation 
value goes to zero. There are three such singular points 
which are related to each other by threefold symmetry, 
and their exact location can be calculated in the k • p 
theory. Note that although the a, diverge in the limit 
k z —> 1 — £, the points where k z = 1 — £ correspond to 
the degenerate band A$, and since ax, ft and a 2l j3 2 have 
opposite signs, the total spin (S)i + (S) 2 = for states 
in A 6 . 

We have carried out ab initio calculations to con- 
firm that tensile strained HgTe and Hgo.5Zno.5Te ex- 
hibit identical spin texture at the Fermi surface. Fig- 
ure [8] illustrates the spin texture on the Fermi surface 
of Hgo.5Zno.5Te. The spins are locked in a two dimen- 
sional plane perpendicular to the (111) axis, and their 



winding number changes from one end of the Fermi el- 
lipsoid to the other. Furthermore, the abrupt change 
in the winding number is accomodated by the vanishing 
of the spin expectation value at three rotation symme- 
try related points on the Fermi ellipsoid. Density func- 
tional theory calculations were performed in the QUAN- 
TUM ESPRESSO suite using the Perdew-Burke-Ernzerhof- 
type generalized gradient approximation (GGA)P^. All 
calculations used an 8 x 8 x 8 Monkhorst-Pack k-point 
mesh with a plane-wave cutoff of 50 Ry. The pseudopo- 
tentials representing the atoms in these simulations were 
generated by the OPIUM package; the pseudopotentials 
were norm-conserving and optimized, and included the 
full relativistic correctio n 1 17 * 18 !. 

We have shown that HgTe tensile strained in the (111) 
axis is semimetallic and exhibits a Fermi surface con- 
sisting of two ellipsoids that contact at a point. Due to 
mirror and threefold symmetry of the zincblende lattice, 
the spin texture on the Fermi ellipsoids is locked in a 
two dimensional plane and vanishes at special singular 
points. We have derived a k • p theory to understand 
this phenomenon at a qualitative level. We propose that 
this symmetry protected spin texture can be observed in 
Hgi_ x Zn x Te for values of x small enough to allow band 
inversion while sufficiently breaking symmetry. 

Note: We recently became aware that the authors of 
Reference UHl studied the effects of biaxial strain (referred 
to as 'strain' in that paper) along the (111) axis in HgTe. 
They found that the fourfold degeneracy at T was lifted 
and no accidental degeneracies resulted for either positive 
and negative strain. Therefore, biaxial strain (compres- 
sive or tensile) is similar to compressive uniaxial strain in 
terms of its effect on band structure. In contrast, we find 
that uniaxial tensile strain not only maintains the acci- 
dental degeneracies already present in HgTe but shifts 
them further away from T. 
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FIG. 6: Spin texture on the Fermi surface of tensile strained HgTe as predicted by the k • p theory of Eq. [T] Top Left: Band 
dispersion in the k • p theory as a function of momentum along the (111) axis. The Fermi level £f (dashed line) is set to be 2 
eV in the units of Eq. [l] Bottom Left: The Fermi surface for £p — 2 eV at the k- p level consists of two spherical shells which 
touch at the point k z = — 1. Top Right: The spin texture of valence states on the blue Fermi spherical shell as a function of 
momentum in the k x , k y plane. Each circle corresponds to a constant k z (indicated by dashed yellow lines cutting through the 
Fermi surface on the left). Bottom Right: The spin texture of valence states on the red Fermi spherical shell as a function of 
momentum in the k x , k y plane. The spins are locked in the k x , k y plane for both Fermi shells. The winding number of spins on 
each Fermi shell changes as one moves from one end along k z to the other. There are three points indicated in yellow where 
the spin texture vanishes. The spin winding number flips from ±2 to =pl at these spin singular points. 
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FIG. 7: (a) Band structure of Hgo.5Zno.5Te is the same as tensile strained HgTe [Fig. 3(b) except the shape and size of 



the Fermi surface and the Fermi velocities in the two materials, (b) Band dispersion at the Fermi level near the accidental 
degeneracies along the (111) axis. There is a twofold degenerate band (thick line) intersecting two nondegenerate bands (thin 
lines), essentially identical to tensile strained HgTe [Fig. |4j. 
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FIG. 8: Spin texture on the Fermi surface of Hgo.5Zno.5Te at the ab initio DFT level, (a) Spins on a Fermi surface ellipsoid of 
Hgo.5Zno.5Te are locked in the plane perpendicular to the (111) axis. The spin winding number around the (111) axis changes 
from +1 to -2 from the left end of the ellipsoid to the right end. (b) Spin texture in the two dimensional plane transverse to 
the (111) direction. Each concentric loop of spins corresponds to a fixed value of momentum along the (111) axis. The spin 
winding number around the (111) axis changes from the inner circles to the outer circles on the Fermi ellipsoid. Since spins are 
locked on the transverse plane, their winding number changes as they vanish at three rotation symmetry related singular points 
(drawn as white circles outlined in black). Left: A higher resolution image showing the change in the spin winding number 
from +1 to -2 as we move from the inner radii to the outer radii on the right. For visual clarity, the magnitude of all vectors 



is normalized in (a) while in (b) the size of each vector equals the square root of its magnitude. 
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